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Abstract 

A non-perturbative method which combines constraints from chiral symmetry break- 
ing and coupled channel unitarity is used to describe the meson-meson interaction up 
to about 1.2 GeV. The approach uses the 0(p 2 ) and 0(p 4 ) chiral Lagrangians. The 
seven free parameters of the 0(p 4 ) Lagrangian are fitted to the data. The results are in 
good agreement with a vast amount of experimental analyses. The amplitudes develop 
poles in the complex plane corresponding to the /o, ao, p, K* , (ft, a and k resonances; the 
latter two, very broad. The total and partial decay widths of the resonances are also 
well reproduced. Further extensions and applications of this chiral non-perturbative 
scheme are also discussed. 
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1 Introduction 



The meson-meson interaction has been the key problem to test Chiral Perturbation Theory 
(xPT), which has proved rather successful at low energies H. The underlying idea is 
that an expansion in powers of the meson momenta converges at sufficiently low energy, 
which in practice is ^/s < 500 MeV. However, the convergence at higher energies becomes 
progressively worse. Even more, one of the peculiar features of the meson-meson interaction 
is the presence of resonances like the a, fo, a® in the scalar sector and the p, K* or the 4> in 
the vector channels. These resonances will show up in the T matrix as poles that cannot 
be obtained using standard X-PT. Nevertheless, the constraints imposed by chiral symmetry 
breaking are rather powerful and not restricted to the region where X-PT is meant to converge 
I- 

Two independent approaches of non perturbative character have extended the use of 
chiral Lagrangians to higher energies and have been rather successful, reproducing important 
features of the meson-meson interaction including several resonances. Although these two 
approaches look in principle rather different, they share a common feature which is the 
imposition of unitarity. One of them 0, g, based upon the Inverse Amplitude Method 
(IAM), first suggested in ||, makes use of the lowest order, 0(p 2 ), as well as the next to 
leading order, 0(p A ), Lagrangians. Elastic unitarity is imposed and thus no mixture of 
channels is allowed. Then, the coefficients of the 0(p 4 ) Lagrangian are fitted to the data. 
The absence of coupled channels has obvious limitations, but in spite of them, the IAM is 
able to generate dynamically the p, K* and a resonances, and to reproduce tctt scattering 
in the (I,J) = (0,0), (1,1), (2,0) partial waves, as well as in the (3/2,0), (1/2,1) and (1/2,0) 
channels of irK scattering. The results are very successful up to 1 GeV in all these channels 
but the (0,0), where it only yields good results up to 700 MeV. The limitations of this single 
channel approach become evident, for instance, in the / (980) and a (890) resonances (J=0 
and 1=0 and 1, respectively) which do not appear as poles in the T matrix. The method 
also has a pathological behavior close to the T matrix zeros 0. 

The second approach dealt with the J=0 sector alone [§]. The input consists of the 0{p 2 ) 
Lagrangian, which is used as the source of a potential between mesons. This potential enters 
in a set of coupled channel Lippmann-Schwinger (LS) equations (actually closer to Bethe- 
Salpeter equations, since relativistic propagators are used) which leads to the scattering 
matrix. The method imposes unitarity in coupled channels; hence it yields inelasticities 
when inelastic channels open up. Amazingly, the approach has only one free parameter, 
which is a cut-off that regularizes the loop integrals of the LS equation. Such a method 
proves rather successful since phase shifts and inelasticities are reproduced accurately up to 
1200 MeV. The /o(980) and ao(980) resonances appear as poles of the T matrix for 1 = 
and 1, respectively, and their widths and partial decay widths are very well reproduced. In 
addition, one finds a pole when I = at y/s ~ 500 MeV with a width of around 400 MeV, 
corresponding to the a meson, which was also found with similar properties with the IAM 

The appearance of the f and a® is due to the introduction of the KK channel in addition 
to 7T7T in I = and irt] in I = 1. These resonances disappear if the KK channel (not considered 
in ||, H) is omitted, while the a in I = is almost not affected. This explains why the f 
and a resonances did not show up in the IAM 0, |5| . 

The success of the scheme of ref. || in the scalar sector gives hopes that it could be used 
in other channels. However, one soon realizes that it does not reproduce properly the J = 
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1 sector. This looks less surprising when one recalls that the 0(p 4 ) chiral Lagrangian can 
be reproduced with the resonance saturation hypothesis ||. That is, assuming that the ac- 
tual values of the 0{p i ) parameters are basically saturated by resonance exchanges between 
Goldstone bosons. In this way, one establishes a clear relation between the information con- 
tained in the 0(p A ) Lagrangian and the resonances in the meson-meson sector, particularly 
vector meson resonances, where the approach of has its stronghold. Indeed, the absence 
of the p and K* in the approach of H], which only uses the 0(p 2 ) Lagrangian, is an indirect 
confirmation of the link between these resonances and the 0(p A ) Lagrangian. 

The approaches of f|, || and || seem complementary and one may wonder whether there 
is a generalization of these methods, containing both them as limiting cases. An affirmative 
answer to this question was recently found and such a generalized method was proposed 
in [|IIJ. The purpose of the present paper is to exploit the idea of [ID] and obtain all the 



predictions of such approach in the meson-meson sector, like phase shifts, inelasticities, 
resonance properties, etc... At the same time we will establish the links between this scheme 
and xPT at low energies. We also illustrate qualitatively, using a toy model, why the 
proposed method is so successful when dealing with amplitudes dominated by resonances. 

2 Unitary amplitude in coupled channels 

Let us write the partial wave decomposition of the meson-meson amplitude with definite 
isospin / as 

T i = Zj(2J + 1)T ij Pj(cos6) (1) 

where T/j is the partial wave amplitude with isospin / and angular momentum J. In each 
one of these channels there are several meson-meson states coupled to each other. In Table 
I, we have listed these states for the J = 0, 1 channels, which contain the most relevant 
meson-meson information below 1 GeV. Note that it is enough to take into account one or 
two states in each channel since we are neglecting here, on the one hand, multipion states 
which are only relevant for higher energies and, on the other hand, the r\r\ that appears for 
(/, J) = (0, 0). The influence of this state is rather small. We have checked it following the 
scheme of || and, although not zero, we found it small enough to omit it with the consequent 
simplicity in the general formalism. 

Hence, throughout the present work, Tjj will be either a 2 x 2 symmetric matrix when 
two states couple, or just a number when there is only one state. In what follows we omit 
the /, J labels and use a matrix formalism, which will be valid for the general case ofnxn 
matrices corresponding to n coupled states. 

The normalization of T is such that 

— 1 k f \ T I 2 (9) 
dtl 64tt 2 s h 1 1/1 1 ' 

where ki and kf are, respectively, the CM three momenta of the initial and final state and s 
is the usual Mandelstam variable. Note that we have chosen a convention for the sign of T 
such that in an elastic amplitude ImT < 0. 
Unitarity in coupled channels implies 

Im T if = T in a nn T* f (3) 
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where a is a real diagonal matrix whose elements account for the phase space of the two 
meson intermediate states n which are physically accessible. With the normalization that we 
have chosen, a is given by the imaginary part of the loop integral of two meson propagators 
in the n state 

<7 nn {s) = lmG nn (s) 
r / \ . f _d^q 1 1 , s 

^nn\S) - I J ( ^ )4 ^ _ ^ + ^ (p _ ^ _ ^ + ^ [4) 

k 

ImG nn (s) = - n / - 6(s - {m Vn + m 2 „) 2 ) 

07T\/S 



where k n is the on-shell CM momentum of the meson in the intermediate state n, P is the 
initial total four-momentum and mi„, r«2 n the masses of the two mesons in the state n. An 
analytical expression for G nn (s) using a cut-off (q m ax) regularization in the integral over d 3 q 
is shown in Appendix A. 

From eq.(|3|) we can extract a and express it, in matrix form, as 

ImG = T _1 • ImT • T* _1 

= -(T- u - T~ l ) = -ImT- 1 (5) 

Hence, 

T-^RcT- 1 -ilmG 
T=[ReT- 1 -ilmG}- 1 (6) 

This is a practical way to write the unitarity requirements of eq.(|3|) which tells us that 
we only need to know ReT -1 since ImT -1 is given by the phase space of the intermediate 
physical states. 

The next point is to realize that the T matrix has poles associated to resonances, which 
implies that the standard perturbative evaluation of xPT will necessarily fail close to these 
poles. As a consequence, one might try to exploit the expansion of T _1 , which will have zeros 
at the poles of T, and in principle does not present convergence problems. For illustrative 
purposes, we can use an analogy with the function tan a; when expanded around x = (x 
playing here the role of p 2 in the chiral expansion). This function has a pole at x = 7r/2. 
Its inverse, cotx, has a Laurent expansion around x = and a zero at x = tt/2. However, 
inverting the expansion of cotx around x = for values of x near tt/2, provides a faster 
convergence than expanding directly tana; around that point. With this idea in mind let us 
expand T _1 in powers of p 2 as one would do for T using xPT: 

T~T 2 + T 4 + ... 

T- 1 ~ T 2 X • [1 + T 4 • T 2 1 ...]- 1 ~ T 2 X ■ [1 - T 4 • T 2 1 ...] (7) 
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This expression requires the inversion T 2 which might not be invertible, as it happens, for 
instance in the (1, 1) channel. In order to avoid the use of T 2 _1 we modify eq.© by formally 
multiplying by T 2 • T 2 -1 on the right and T 2 l • T 2 on the left. All the steps are justified using 
the continuity of the functions involved in the derivation, starting from a matrix close to T 2 , 
which can be inverted. Thus, eq.(^]) can be rewritten as 

T = T 2 • [T 2 • ReT" 1 • T 2 - iT 2 • Im G ■ T 2 }- 1 ■ T 2 (8) 
Now, using the expansion for T _1 of eq.(^|) we find 

T 2 -ReT- 1 -T 2 ~T 2 -ReT 4 + ... (9) 

and recalling that 

ImT 4 = T 2 -ImG-T 2 (10) 
we finally obtain, within the 0(p 4 ) approximation 

T = T 2 -[T 2 -T 4 }~ 1 -T 2 (11) 

Note, as it is clear from eq.(H), that what we are expanding is actually T 2 • ReT -1 • T 2 , 
which in our analogy would be equivalent to x 2 cot x, which is also convergent around x = 0. 

In another context, the above equation can also be derived using Pade approximants [|TTJ. 
This equation is a generalization to multiple coupled channels of the IAM of ref.[|], ||. It 
makes the method more general and powerful and also allows to evaluate transition cross 
sections as well as inelasticities. 

The coupled channel result has additional virtues with respect to the single channel IAM. 
Indeed, in this latter case the expansion of eq.(|7|) is meaningless if |T 2 | < |T 4 | or T 2 = 0. 
In particular, if T 2 vanishes, eq.([TT|) yields T = T 2 ■ T 4 _1 , which has a double zero, whereas 
the correct result would be T ~ T 4 . This indeed occurs in the J = partial waves below 
threshold (Adler zeros). However, within the coupled channel formalism, if a matrix element, 
say (T2)n, vanishes, it is sufficient that (T 2 ) 12 7^ 0, since then eq.(|TTD gives (T) n ~ (^4)11, 
which is the correct result. In conclusion, while the single channel IAM gives a double zero 
whenever T 2 = 0, the coupled channel method leads to single zeros close to the zeros of T 2 . 

The single channel IAM has another related problem, since close to the Adler zero it 
presents an spurious pole when T 2 = T 4 . The coupled channel method also avoids this 
problem, although it runs into a similar one when the determinant of the T 2 — T 4 matrix 
vanishes below threshold. This happens indeed for J = 0, / = around y/s ~ 120 MeV. 
Excluding the neighborhood of this zero of the determinant, we can still recover from eq. (|TT|) 
the usual xPT expansion, T ~ T 2 + T 4 + ... valid for low energies, typically \yfs\ < 500 
MeV. In any case we concentrate here on results above the two pion threshold. 



It is now important to realize that eq. (|lTD requires the complete evaluation of T 4 , which is 
rather involved when dealing with many channels, as it is the case here. Instead, we present 
a further approximation to eq.flTTD which turns out to be technically much simpler and rather 
accurate. In order to illustrate the steps leading to our final formula, let us make before 
another approximation. Let us assume that through a suitable cut-off we can approximate 

Re T 4 ~ T 2 • Re G ■ T 2 (12) 
In such a case we go back to eqs. (]8|, ||) and immediately write 
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T=[l-T 2 -G\- 1 -T 2 (13) 

that is equivalent to 

T = T 2 + T 2 -G-T (14) 

which is a LS equation for the T matrix, where T 2 plays the role of the potential. This is 
actually the approach followed in ref. ||. 

There is a subtle difference between eq.(O) and the ordinary LS integral equation. In- 
deed, eq . fli~4]) is an algebraic equation since T 2 and T are factorized out of the integrals with 
their on-shell value. In contrast, in the ordinary LS equations, the T 2 GT term is actually 
the integral of eq.(|4]), including T 2 and T inside the integral, since both of them depend on 
q. Due to the structure of the 0{p 2 ) Lagrangian, it was shown in || that writing T 2 (q) as 
T 2 ons/ie "(g) + T 2 ^ shell (q), the off-shell part renormalizes couplings and masses and hence it 
had to be omitted. Therefore T 2 , and T factorized outside the integral with their on-shell 
values. As a consequence, the very same algebraic equation (|TJ) was obtained. 

As we have already commented, the approximation of eq.(0) leads to excellent results 
in the scalar channels. However, as we mentioned in the introduction, the generalization to 
J 7^ is not possible since basic information contained in the 0(p A ) chiral Lagrangian is 
missing in eq.([T^). The obvious solution is to add a term to eq.(|T2p such that 

ReT 4 ~T 4 p + T 2 -ReG-T 2 (15) 

where T 4 P is the polynomial tree level contribution coming from the 0(p A ) Lagrangian, whose 
terms contain several free parameters, usually denoted Li. Within our approach, these 
coefficients will be fitted to data and denoted by Lj since they do not have to coincide with 
those used in xPT , as we shall see. Actually, the Li coefficients depend on a regularization 
scale (/i). In our scheme this scale dependence appears through the cut-off. 

In addition, there are also differences between our renormalization scheme and that of 
standard xPT . Indeed, our approach considers the iteration of loop diagrams in the s- 
channel, but neglects loops in the u or t channels. However, the smooth structure of these 
terms for the physical s-channel, since we are far away from the associated singularities, 
allows them to be approximately reabsorbed when fitting the L\ coefficients. Concerning 
tadpoles, they would be exactly reabsorbed in the Li in the equal mass case. Therefore, 
when masses are different, we are omitting terms proportional to differences between the 
actual masses squared and an average mass squared. Thus all these contributions will make 
the Li differ from the L i; although we expect them to be of the same order. 

This way of dealing with tadpoles has an additional advantage. Apart from the usual 
tadpole diagrams that would also appear in standard xPT there are some additional tadpole 
terms. They come from the determinant of the SU(3) metric that should be included in the 



path integral measure in order to make the generating functional SU(3) covariant [I3"|. With 
dimensional regularization such contributions vanish, but that is not the case when using a 
cutoff regularization [TJ|]. Nevertheless, we have just described how tadpoles are absorbed 



within our approximation and thus we do not have to calculate them. 

With these approximations our calculations have been considerably simplified at the ex- 
pense of losing some precision at low energies with respect to the full 0(p A ) xPT calculation. 
As far as we are mostly interested in resonance behavior as well as higher energies this is 
not very relevant. Nevertheless, if the complete 0{p A ) calculations were available, we could 
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directly use eg. (|TT|) , and have both an accurate low energy description and a good coupled 
channel unitarity behavior. 

Using eqs.(|§|) and (Pf), our final formula for the T matrix is given by 



T = T 2 -[T 2 -T[ -T 2 -G- Ta]- 1 • T 2 



(16) 



3 Toy model 

In order to illustrate how the method works, we take a simple case of one channel and one 
amplitude around a resonance which we assume to know exactly. That is, 



ap 



(17) 

q 2 — M 2 + i 2MV V ; 

where p 2 is an invariant quantity, of dimension momentum squared, related to the momenta 
or masses of the pseudoscalar mesons, q the total four-momentum of the meson pair and 
2MT = —ap 2 lmG. The above equation satisfies unitarity exactly as can be seen by using 
eq.(|). 

To Oik 2 ), k = p, q, we have 



whereas at 0{k A ) we have 



ReT 4 



Then, using eq. (|TTD we find 



2 2 

ap q 

' M 4 



' M 2 



(19) 



T-- 



rp2 

1 2 



ap 



T 2 — Re T 4 — iT 2 Im GT 2 
ap 2 

q 2 — M 2 — iap 2 lm G 



M 2 (l 



i 

M 2 



ia 



Af 



■ImG) 



(20) 



So, as we can see, in this particular case the IAM leads to the exact result for the T 
matrix, eq. (|l7D . The result is exact here because T 2 ■ ReT -1 ■ T 2 is an 0(k 4 ) function and 
hence the expansion up to 0{k A ) in eq.(|]) is exact. However, the structure of eq. (|T7|) is that 
of a meson propagator of an unstable particle like the fo,a , p, K* , etc... resonances. This 
could justify why the scheme which we propose works even better that one could naively 
anticipate, at least for resonant channels. 

The above argumentation uses the same power counting in momenta as x^T, but pre- 
sumes that the 0(k 2 ) amplitude comes from the exchange of a resonance. This seems to 
be in conflict with 0, where it is shown that resonance exchange contribution shows up at 
higher orders. However, when taking into account requirements of short distance behavior 
of QCD, these two points can be reconciled. In fact, this has been shown, in [I5|, where 
a classical vector meson dominance expression for the pion form factor is obtained, in the 
same lines as eq.(|l7]), starting from chiral Lagrangians and imposing those QCD constraints 
at short distances and the large iV c limit. 
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In relation to the previous arguments, the link between unitarized xPT and vector meson 
dominance has also been discussed in fill . 



4 The matrix elements of To and Tk. 



The lowest order chiral Lagrangian is given by 



f 2 

C 2 = L.(d^U j d»U + M(U + U ] )) 



(21) 



where / is the pion decay coupling and () stands for the trace of the 3x3 matrices build 
out of £/($) and M. 



U($) = exp(z^2$//) 
where $ can be expressed in terms of the Goldstone boson fields as 



7T 



K + \ 



V 

The mass matrix M is given by 



7T" -^7T° + ^ K° 

K- 



K° 



75V 



M 



(ml 
ml 

2 ^2 



\ 2m K — ml 

where we have assumed the isospin limit m u = m^. 
The 0(p 4 ) Lagrangian is given by 



(22) 



(23) 



(24) 



C^Lxid^d^U) 2 + L 2 {d^d v u)(d^dru) 

+ Lz(d^U ] d»Ud u U ] d v U) + U(d^U ] d»U)(U ] M + M*U) 

+ L 5 (d^d^U(U + M + M + U)) + L 6 (U^M + M + U) 2 

+ L 7 (U ] M - M ] U) 2 + L 8 (M^UM^U + U ] ' MU ] M) (25) 

where the terms which couple to external sources are omitted [|TJ. 

The states with definite isospin, with the phases \tt + ) = — 11,1), \K~) = — |l/2 — 1/2), 
are given by 

1 = 0, 

\KK) = -l=\K + {q)K-{-q)+K\q)K\-q)) 

i = i,i 3 = o, 

\KK) = -^\K + {q)K-{-q)+K\q)K\-q)) 
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7T 



rj) = \ir°(q)r)(-q)) 



Itttt) = --| 7r +(g) 7r (-g)-vr (q)n + (-q)) 
I = 2,I 3 = 2, 

k^) = ^k + (<?> + (-<?)) 

I = l/2,I 3 = l/2, 

^ _ 1 



\Kir) = 
\Krf) = \K+(<f)r](-<f)) 



7r + (g)ir (-g) + -i= 7 r (g)^ 



-5)> 



J = 3/2, J 3 = 3/2, 

|^tt) = -|^ + (g> + (-g)> 

We should note that in the states of identical particles we have included an extra 1 / \[2 
factor in the normalization. This is done to ensure that the resolution of the identity gives 
unity (recall that E ? |7r°(g)7r°(— g))(7r°(g)7r°(— g)| = 2 with the states 7r°(g)7r°(— q) normalized 
to unity). This normalization yields the ordinary unitarity formulae, eq.(|3]), which we are 
using to extract phase shifts and inelasticities. However, we should return to the proper 
normalization at the end in order to obtain the physical amplitudes. 

The amplitudes which we obtain are compiled in Appendix B. The projection over each 
partial wave J is done by means of 



T, 



i j 



P J (cos9)T I (6)d(cost 



(26) 



In the case of two coupled channels, T/j is a 2 x 2 matrix whose elements, (Tjj)ij are 
related to S matrix elements through the equations (omitting the /, J labels) 



CO n 

COia 



-^I(5)u-1], (T) M 
(T) 21 - -Ja^L 



2ip 2 



[(5)22 " 1] 



2i v /pip2 



(S) 



(27) 
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with pi, p2 the CM momenta of the mesons in state 1 or 2 respectively. The S matrix has 



the structure [17 



S 



i(l-rf) 



rje 2iSl 

2M/2, 



2M/2 e i{&i 



D i(5i 



t(l-T] 2 ) 

-8 2 ) rn P 2i8 2 



S 2 ) 



(28) 



rje 

where 8\ and 82 are the phase shifts for the elastic 1 — > 1 and 2^2 processes (for instance, 
KK —>■ KK and mi — > Tin in (J, J) = (0, 0)) and 77 is the inelasticity. 

It is interesting to note that, by means of (T)n and (T) 2 2 one can determine t],8i and 
82, and hence the (T)i2 = (T) 21 matrix elements are redundant. We determine them from 
our coupled equations and verify that the structure of eq.(p7|) is satisfied, which is another 
check of the coupled channel unitary that we have imposed from the beginning. 
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5 Results 



We have carried out a fit to the data, which is shown in figs. 1 to 7, using as free parameters 
the Li with i = 1,2,3,4,5,7 and 2L 6 + L 8 . The cut-off is fixed to q max = 1.02 GeV. The 
values which we obtain are shown in Table II. By comparing them with the standard values 
for the Li coefficients obtained in xPT at the scale p = 2q max /^/e (see appendix A. 2) we 
see that they are of the same order. 

We show first the results on phase shifts and inelasticities in the different channels and 
later on we discuss about the pole positions, widths and partial decay widths. 

5.1 Phase shifts and inelasticities 

We will now go in detail through the results in each (I, J) channel. 

5.1.1 Channel (0,0) 

As we can see in eq.(|27D we have three independent magnitudes Si, 5 2 and r\. In figs. (1. a) and 
(l.c) we show the b\ and 5 2 corresponding to tttt — > tttt and KK — > KK elastic scattering. 
In fig.(l.b) we plot the phase shift for KK — > tttt. This is actually 5i + 5 2 , which is therefore 
redundant information. However, there are data for this process but not for elastic KK, 
and that is why we are plotting 5i + 5 2 . The agreement with experiment is good, with small 
discrepancies in the KK — > tttt phase shifts. In fig.(l.a) we see a bump around 600 MeV 
which is due to the a resonance, whose associated pole appears around 442 — i225 MeV, as 
we shall see below. The fast raise in the phase shift at 1 GeV is caused by the fo pole around 
980 — H4 MeV, which translates in an apparent mass of ~ 980 MeV and a 30 MeV width. 
Small discrepancies with data start showing up around 1.2 GeV. The omission of the rjr] and 
four meson states should limit the validity of the approach at high energies since then these 
channels start being relevant. 

5.1.2 Channel (1,1) 

In fig. (2. a) we display the tttt — > tttt phase shifts which clearly show the p meson. The perfect 
coincidence of the results with the very precise data indicate that both the position and the 
width of the p are very well described. In fig.(2.c) we show the phase shifts for KK — > KK 
scattering, for which there are no data. As we can see, they are very small, which implies a 
weak KK interaction. Therefore the 5\ + 5 2 phase shift of KK — > tttt is essentially that of 
tttt — > tttt. The fact that the inelasticity is practically one, indicates that there is almost no 
mixture of tttt and KK. This feature makes the p to behave as a pure tttt elastic resonance. 
That is why the single channel IAM gave essentially the same results as obtained here [|J. 

5.1.3 Channel (2,0) 

The I = 1 tttt scattering contains only one state as shown in Table I. In fig. (3) we show the 
resulting phase shifts, whose agreement with experimental data is remarkably good up to 

1.2 GeV 
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5.1.4 Channel (1,0) 



In fig. (4. a) the irt] — > 7r^ phase shifts are shown. Those of KK — > KK are plotted in 
fig.(4.b) and the inelasticities in fig.(4.c). In the latter, it can be seen that there is an 
appreciable mixture between ttt] and KK above KK threshold. In fig.(4.d) we compare 
a mass distribution for tti] around the region of the ao resonance. The data are obtained 
from |34J] using the K~p — > S + (1385)7r^ reaction, whose cross section (following [j35] ) can be 



written as 

£ - c ^ < 29 > 

where m is the 7T~r] invariant mass, q the tt momentum in the n^i] CM frame, t the 
n^i] — ► tt~i] scattering amplitude and C a normalization constant. We observe a fairly 
good agreement with the experimental numbers. 



5.1.5 Channel (1/2,0) 

The two coupled states are now Kir and Krj. In Fig. (5. a) we plot the phase shifts for 
Ktt — > Kir. The agreement of the results with the data is quite good, although a bit on the 
upper part. The results and the data show a broad bump, which is related to the presence 
of a pole which appears around 770 — 2250 MeV. Such a resonance, whose existence has been 
claimed in a recent data analysis [O], is predicted in quark models of q 2 q 2 systems K| and 



is usually denoted by /t(900). This resonance bears some similarity with the a in the (0,0) 
tttt elastic scattering channel, which is also very broad. Finally, the Kr] — > Kr] phase shifts 
are small as shown in fig.(5.c) and the inelasticities given in fig.(5.d) are not distant from 
unity. This fact indicates a small mixture of Ktt with Krj. 



5.1.6 Channel (1/2,1) 

In this case we also find a resonance in fig. (6. a), analogous to the p, but in the Ktt system. 
This resonant state, known as the if* (892), is as clean as the p, and the agreement of our 
results with the data is remarkably good over the whole range of energies up to 1.2 GeV. In 
fig.(6.c) we plot the Kr] — > Krj phase shifts, which are very small. Finally, in fig.(6.d) we 
can notice that rj « 1 which means that there is practically no mixture of Kir and Krj in 
this channel. This justifies the success of []4j reproducing this resonance using only the Ktt 
state and elastic unitarity. 



5.1.7 Channel (3/2,0) 

In fig. (7) we show the Ktt phase shifts. As we can see in the figure, the agreement with the 
data is quite good up to about 1.2 GeV. 

The channel (3/2,1) in Ktt (see table I) is such that T 2 = 0, since there is only S-wave 
there. In this case our method cannot be applied, as discussed above, and we should just 
take the T 4 contribution. That also happens for the J = 2 channels, since the structure of 
T 2 , which is 0(p 2 ), is a linear combination of s, t, u and squared masses. Therefore there is 
only J = 0, 1 in T 2 , but not J = 2. Hence, the lowest contribution can only be obtained from 
the T4 terms and our method has nothing to improve there with respect to xPT . The phase 
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shifts in these channels are small and have been discussed in Hence we omit any further 
discussion, simply mentioning that the agreement with data found in |4j is fairly good. 

There is another interesting result in the (0,1), channel which is the appearance of a pole 
around 990 MeV, that we show in fig. (8). Below 1.2 GeV there are two resonances with such 
quantum numbers. They are the u and the 0, which fit well within the qq scheme, with 
practically ideal mixing, as -^(uu + dd) and ss, respectively. Hence, the uj would almost 
decouple from KK and then we should not expect it to appear in our scheme with only the 
KK channel. The three pion channel, into which the uo mostly decays, is not considered in 
our approach, restricted to two meson states. In contrast, the couples strongly to the KK 
system. It seems then natural to identify the above mentioned pole with the 0. Indeed, it 
is only 30 MeV below its real mass, 1020 MeV (which means a relative deviation of only 3 
%). Due to this shift towards lower energy this resonance appears below the KK threshold 
and this is why we find no width at all. Nevertheless, as far as its physical width is only 
~ 4 MeV, it seems plausible that that the small coupling to three pions (an OZI suppressed 
coupling of third class) which we are not taking into account, could be enough to improve 
the agreement between the position of our resonance and its real mass and width. 



5.2 Pole positions, widths and partial decay widths. 

We will now look for the poles of the T matrix in the complex plane, that should appear in 
the unphysical Riemann sheets ( the conventions taken are those of ||, which can be easily 
induced from the analytical expressions of Appendix A). Let us remember that the mass 
and the width of a Breit-Wigner resonance are related to the position of its complex pole 
by y/s po i e ~ M — iT/2, but this formula does not hold for other kind of resonances. In table 
II we give the results for the pole positions as well as the apparent or "effective" masses 
and widths that can be estimated from phase shifts and mass distributions in scattering 
processes. Note that such "effective" masses and widths depend on the physical process. 

We shall make differentiation between the p and K*, which are clean elastic Breit-Wigner 
resonances, and the rest. For the p and K* their mass is given by the energy at which 5 = 90° 
and the width is taken from the phase shifts slope around 5 = 90°, by means of 

r ^ = ^r tan5(s) (30) 

We also saw that, in practice, the p and K* only couple to tttt and Kit, respectively. The 
a decays only to tttt and the k only to Kit due to phase space and dynamical suppression 
of other channels (see fig.(5d)). The case of the f and a is different, since they can decay 
either to tttt or KK (the fo), and ttt) or KK (the ao). In order to determine the partial decay 
widths of these resonances we follow the procedure of f|, where we show that, assuming a 
Breit Wigner shape for the amplitudes around the resonance pole, the partial decay widths 
are given by 



j r ^max Q 

1D7T JE min Hi 



1 f E max q (ImT 21 ) 2 , , 

167r 2 JE min E 2 Im Txi 

where E stands for the total CM energy of the meson-meson system, q is the momentum of 
one meson in the CM and the labels 1, 2 stand for KK, tttt in the case of the fo and KK, 



12 



7r?7 in the case of the do- The masses of the final mesons are mi, m 2 . The upper limit in 
the integral, E max , is ~ M R + T R where T R is the total width f8[ and E min = M R — T R , 
unless the threshold energy (mi + m 2 ) for the decay is bigger than that quantity, in which 
case E min = mi + m 2 . In this way we largely avoid the contribution of the backgrounds in 
the amplitudes. One caveat must be raised concerning eq.([H]), which was already pointed 
out in the study of the f — > 77 decay [[43|]. The subtlety is that around this resonance the 
phase shifts (see fig.(l.a)) are of the order of 90°, due to the background coming from the 
broad a pole. This background makes the fo — > nn coupling constant to appear effectively 
multiplied by a tt/2 phase (i factor) and in this way the T12 amplitude around the fo looks 
like an ordinary Breit Wigner multiplied by i. This means that the real part has a peak 
around the resonance and the imaginary part changes sign. In this case the arguments used 



in H and [|3 lead to a trivial modification in T R2 , where ImT 12 should be substituted by 
ReT 12 . 

It is also very instructing to see the representation of the poles in a three dimensional 
plot. In fig. (9) we are showing on the left the imaginary part of the (0,0) txti — > tttt scattering 
amplitude on the second Riemann sheet. It is possible to see very clearly the appearance 
of two poles that correspond to the a and the fo resonances. The former is located at 
442 — i227 and thus is very far away from the real axis, which implies a huge effective width. 
In contrast, the other pole is located at 994 — il4 MeV accordingly to the narrow width of 
the fo resonance. 

Apart from the position of the poles, there is an additional piece of information which 
also determines the observed shape of a resonance. It also explains some of the differences 
between the "effective" masses and the real part of the pole position. On the right of fig. (9) 
we give a contour plot, again of the imaginary part of the (0,0) amplitude in the second 
Riemann sheet. Notice that both poles are oriented differently, almost transversally, on the 
complex plane. On the one hand, the fo pole is oriented almost perpendicularly to the real 
axis, which is the relevant one in this work. As a consequence, in the positive real axis, the 
imaginary part of the amplitude first grows rapidly and then drops very fast again, giving 
rise to the dramatic variation of the phase shift typical of resonances. A similar orientation 
is found for the p, K* and a>o resonances too. On the other hand, the a pole is oriented so 
that in the real axis we only see a slow and smooth increase, but almost no decrease, of the 
imaginary part. That is also the case of the k resonance. This feature, together with the 
fact that both the a and the k are very far from the real axis explains why it is so hard to 
establish firmly their existence and their physical parameters. 

Finally, in fig. (10) we present a very detailed contour plot of the p and ao poles. Both of 
them are almost perpendicular to the real axis, but the former is tilted clockwise, whereas the 
latter is tilted anti-clockwise. Let us now remember that the real part of the pole position, 
roughly, should give us the apparent mass of the resonance. However, the lines of maximum 
gradient of each pole cross the real axis at a point which is slightly different from the real 
part of its position. Therefore, those poles rotated clockwise, as the p or the K*, have an 
apparent mass a little bit higher than that given by the pole position. In contrast those 
tilted anti-clockwise, yield a resonance whose mass is somewhat lower that the one obtained 
from the pole. That is the case of the fo and the a . 
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6 Conclusions and Outlook 



We have used a coupled channel unitary approach, together with the dynamical information 
contained in the 0(p 2 ) and 0(p A ) chiral Lagrangian, which allows us to study the meson- 
meson interaction up to about 1.2 GeV. This non-perturbative method generates poles in 
the complex plane corresponding to physical resonances. We have used the experimental 
information available to make a fit of the 0(p A ) Lagrangian coefficients. These are Li, 
i = 1,2, 3, 4, 5, 7 and 2L 6 + L 8 , whose actual values depend on the cut-off that we have used 
to regularize divergent one loop integrals. With those seven degrees of freedom we are able 
to fit, up to 1.2 GeV, all the experimental information in seven meson-meson channels. Each 
one of this channels consists of two phase shifts and an inelasticity. Moreover, in our results, 
we obtain the position and widths, partial decay widths, etc.... of all the resonances that 
appear in those channels below 1.2 GeV. Apart from the standard f , a , p, K* resonances, 
we find poles in the T matrix for the a in the tttt I = J = channel and for k in the (1/2,0) 
channel, both them very broad. 

The method has proved very efficient to extend the ideas of chiral symmetry at energies 
beyond the realm of applicability of xPT. However, at energies higher than 1.2 GeV, the 
limitations of the model show up, since, among other things, we have restricted ourselves 
only to two meson states. The restrictions in the space of states precluded the appearance of 
the uj resonance which couples dominantly to three pions. However, the resonance which 
couples strongly to KK does appear in the scheme, although slightly shifted towards lower 
energies. Presumably, by including the coupling to three pions, although very small, it 
should be enough to shift the mass to its correct place. 

One of the formal weakness of the approach is that loops in crossed channels, as well 
as some tadpole contributions, are not explicitly included in the calculation. In practice, 
their effect can be reabsorbed in the fit of the 0(p A ) parameters, whose values can then be 
different from those obtained for the standard low energy xPT approach. 

This approximation could be improved by using eq . (|TT|) with the full 0(p 4 ) xPT calcu- 
lation, which includes one loop in crossed channels and the tadpoles. This would allow a 
more straightforward comparison with standard xPT as well as a better accuracy in the low 
energy results. Although such calculations are welcome and there is indeed some work in 
progress j44j, they are far more involved to calculate and use. 



Applications of the method to other physical problems are also in order. Indeed, it can be 
easily extended to deal with processes where meson pairs appear in the initial or final state, 
like meson pair photoproduction ||43|| . It looks likely that it could also prove useful describing 
the meson-nucleon interaction [^5[ complemented with Heavy Baryon Chiral Perturbation 
Theory. In addition, the method, non perturbative in nature, is equally well suited to study 
the meson-meson interaction in a nuclear medium where there has been some speculation 
about the appearance of bound im pairs |46[] . 

Finally it seems that the approach could be extended to the effective chiral Lagrangian 
description of the Standard Model Strongly Interacting Symmetry Breaking Sector, where 
the single channel approach has already been applied . 
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A Analytical formula for G(s): Relation between cut- 
off and dimensional renormalization 

In this appendix we are showing the relationship between our regularization scheme and 
dimensional regularization, which is the usual one when dealing with xPT. 

A.l Analytical formula for G(s) with a cut-off regularization 

In the general case with different masses, Mi and M 2 



G(s) 



32tt 2 



log 



Ml 
Ml 



A 

+2- loj 

s 



M? M.j 
+ log^^ 

Qmax 



+ 



V 



log 



A + v,/l + #- 



s + A + uJl + p- 



-S + A + 1/J1 + 



+ log 



A + 1/4/I 



Ml 



Ml 



2 log 



\ 



Ml 

Qmax 



1 + 



\ 



1 + 



Qmax 



(32) 



where v = ^J(s - (Mi + M 2 ) 2 )(s - (Mi - M 2 ) 2 ) and A = Mf — M|. In the case of equal 
masses, Mi = M 2 = m, the above formula reduces to 



G(s) 



(An) 



alog- 



aJl + 



+ 1 



<tJ1 + 



- 1 



-2 log 



m 



\ 



1 + 



Qmax 



(33) 



where now, a = Jl — Am? Is. 

The numerical evaluation of the principal part of eq.(Q) is also performed as an additional 
check. 



A. 2 Relation between the cut-off and the dimensional regulariza- 
tion scale 

In order to obtain the relationship between the cut-off and the renormalization scale fi 
let us consider, for the sake of simplicity, the case with equal masses (the same result is 
obtained with different masses but the formulas are more cumbersome). As far as we are 
going to compare the same function calculated in different ways, let us denote by G (s) the 
G(s) calculated with a cut-off regularization and G D (s) the one calculated with dimensional 
regularization. In this latter case we have 
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G D (s) 



(47T) 1 



1 , 2 ,£7 + 1 

2 + log m + a log 

e a — 1 



(34) 



where 1/e = 1/e - log(47r) + 7 with D = 4 + 2e. 

The scale \i in G D (s) appears through the inclusion of the at (9(p 4 ) 



A = + TjA (35) 

where L^(fi) is the renormalized value of Li at the /i scale, I\ is just a number and 



A 



32tt 2 



~ ~\~ lOg fi 1 



(36) 



The log/z 2 , and its companion - — 1, are incorporated in G D (s) so that at the end one 
has a logarithm of the dimensionless quantity m 2 //i 2 . In this way we rewrite G D (s) as: 



G D (s) 



(4tt) 



, m z , a + 1 
-1 + log — + a log - 

jLT cr — 1 



(37) 



We expand eq.(|33|) in powers of m 2 /q 2 nax to compare with the cut-off regularization, as 
follows 



G c (s) 



(An) 2 
1 

(4tt) 2 
1 

(i^j 2 




-1 +log 



(38) 



Then comparing eqs. (|37D and one has: 



2g r , 



1.2 q ri 



(39) 



Hence, to our cut-off q max ~ 1 GeV would correspond a /i = 1.2 GeV dimensional 
regularization scale. In Table II, we have listed the values of the L t parameters and those 
of standard xPT scaled to fi = 1.2 GeV. As it is explained in the text, in our fit we have 
neglected the crossed channel diagrams and we have treated tadpoles differently. The effect 
of these contributions is effectively reabsorbed in our Lj parameters, hence some differences 
between the Li and Lj parameters should be expected and this is indeed the case as can be 
seen in Table II. Note that, even if we had used the complete 0(p 4 ) xPT calculations, these 
parameters could be different, since they have been obtained from a fit over a much wider 
range of energies than it is used in xPT and higher order contributions have been included. 

Finally, note that the terms 0{m 2 j q max ) in eq.(|38f) yield 0(p 6 ), or higher, contributions 
and that is why they are not included in G D (s). 

It is also worth stressing that the relationship of eq.(^) is independent of the physical 
process and channel since the function G(s) appears in all them in the same way. 
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B Amplitudes 



We have used the following formulae in our calculations. Note that, as it has been explained 
in the text, we have an overall sign of difference with the definitions in 0, as well as a 1/2 
factor for those amplitudes with identical particles. 

B.0.1 Masses and decay constants 



fir — fo 



fK — fo 



fr) — fo 



1 + ^(L 5 + L 4 ) 
Jo 

Am 



+ 



5m 



K 



1 + ^f(L 5 + 2L 4 



1 + 



Po 
Am 2 

Jo 



fo 2 
Am 2 

JO 



8m 2 K + Ami 



Po 



(40) 



m^ = m, 



7T 



5m; 



1 + — f(2L 6 + L 8 -L 4 -L 5 ) 
Jo 



lQm 2 K 



fo 



(2L e 



La 



m 



K ' 



-m. 



o K 



1 + 



16m 



K 



fo 2 



(2L 6 + L 8 — L 4 — — L 5 



5m; 



/o 2 



(2L 6 - L 4 ) 



where the subscript refers to bare quantities. 



(41) 



B.0.2 mi — > 7T7T scattering 

The definite isospin amplitudes T^> are obtained from just one amplitude T: 

T {0 \s,t,u) = (3T(s,t,u)+T(t,s,u)+T(u,t,s))/2 
T (1) (s, t, u) = (T(t, s, u) - T{u, t, s))/2 
T< a > (s, t, u) = (T(t, s, u) + T{u, t, s))/2 

where T = T 2 + T 4 is given by: 
„ ml — s 



* f2 

J 7T 

T 4 = — ^ {(2Lx + L 3 )(s - 2ml) 2 + U [(t - 2ml) 2 + (« - 2m') 2 
+2(2L 4 + L 5 )m 2 (s - 2ml) + 4(2L 6 + L 8 )m*} 
which have been obtained at tree level from £ 2 and £ 4 , respectively. 



(42) 



(43) 
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B.0.3 Kn — > Kn scattering 

Using crossing symmetry, we can write the 1 = 1/2 amplitude in terms of that with / = 3/2, 



as 



where 



(a, t, u) = (u, t, s) - V/ 2 ) (s, t, u) 



,(3/2) _ s - {ml + m\) 



(44) 
(45) 



T (3/2) = _ 



J IT J j 



K 



{(4Li + L 3 )(t - 2m 2 )(t - 2m 2 K ) + 2L 2 (m 2 + m\ - sf 



+ (2L 2 + L 3 )(m 2 n + m 2 K - u) 2 + 4L 4 (m 2 + mj c )t - ±m\m\ 



+ L 5 



(ml + m 2 K )(ml + m|- - s) - 4m 2 m|- + 8m 2 r m| r (2L 6 + L 8 )| 



once more, they have been obtained, respectively, from £ 2 and £ 4 at tree level. 
B.0.4 KK -> fTZ scattering 

The definite isospin amplitudes can be written just in terms of two: 

T (0) (s, t, m) = T+-+-(s, t, u) + T 5o+ -(s, t, u) (46) 
T« (s, t, u) =T + - + - (s, t, u) - T 5o+ ~ (s, t, u) 

where T + ~ + ~ is the amplitude for K + K~ — > K + K~ , whose respective C(p 2 ) and 0(p 4 ) 
contributions are 

u — 2m 2 K 



T +- 



(s,t,u)-- 



(47) 



Ti + ~ + ~ (s, t, u) = {2L 2 (u - 2m^) 2 + (2L X + L 2 + L 3 ) [(s - 2m^) 2 + (t - 2m 2 K f 

-2um 2 K {2U + L 5 ) + 8m 4 K (2L 6 + L 8 )} 
whereas T 5o+ ~ is the amplitude for K°K° — > K + K~ , which is given by 

t 2 50+ -M,.)=^^ 



(48) 



T 4 50+ -( S , *,«) = - 7r {(4Li + L 3 ){s - 2m 2 K ) 2 + (2L 2 + L 3 )(t - 2m 2 K f + 2L 2 (u - 2m 2 K ) 2 

JK 

2 1 
+2sm^(4L 4 + L 5 ) + -tm^(2L 4 + L 5 ) - 8m 4 K (2L 4 + L 5 - 2L 6 - L 8 ) j 

B.0.5 nix — > fTfT scattering 

Again, we can use crossing symmetry to obtain, from Kn — > Kn, the definite isospin ampli- 
tudes Tf of this process: 



T (°) = ^ (T^ (u, s, t) + (t, s, uj) 

T (D = ^ ( T (3/2) (M)S?t) _ T (3/2) (t;S; 



(49) 
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B.0.6 Krj — > Krj scattering 

This process is pure I = 1/2. We obtain the following contributions to the amplitude: 

, 6ml + 2ml ~ 9t 
T 2 (s,t,u) = - ' 



12 j ' v fK 



T 4 (s,t,u) 



{2(t - 2m 2 K )(t - 2m 2 r) )(12L 1 + 5L 3 ) + [(u -m\- m\f 



(50) 
(51) 



+ (s -m\- m 2 K f \ (12L 2 + L 3 ) + 2(t - 2m 2 K ) [llm 2 K (2L 4 + L 5 ) - m^(L 4 + 3L 5 
+ Am A K [3(2L 4 + L 5 ) + 32(L 6 + L 7 + L 8 )] + 2(t - 2mJ) [9m^(2L 4 + L 5 ) + (3L 4 - L 5 
+ 4m*(16L 7 + 8L 8 - L 5 ) + 6L 5 mJ - i [m^(24L 4 + 7L 5 ) + 2m2(6L 4 - L 5 ) + 9L 5 mJ 
+ 6m;m|(4L 4 + L 5 ) + 2m^(6L 4 + L 5 ) + 2 m^m^ [6L 4 + L 5 - 8(2L 6 + 7L 8 + 12L 7 )]} 

B.0.7 7f?y — > 7T7r scattering 

The 7 = 1/2 amplitude can be obtained as follows: 



T (1/2) (s,t,u) = ^T^ K _ n+ (s,t,u) 



(52) 



The 0{p 2 ) and C(p 4 ) contributions to K°i] K n + are 



T 2 (s,t,M): 

T 4 { S ,t,U)-. 



V6 \8rnt + 3ml + ml - 9t 



3Gfi<f v 



(53) 



2/3 



[3L 3 2(t -ml- m 2 )(t - 2m 2 K ) - (s - m 2 K - m 2 n )(s - m 2 K - m 2 ) 



Vlf 2 

(u - m 2 K - ml)(u -m 2 K - m 2 ) 



(t + m% — m 2 )(7m 2 K - hml) 



+ 4m^(3t - 3ml - m 2 ) + 2(t - 2m 2 K ){m 2 K + ml) + A(m 4 n - m 4 K ) 
+ + 16(2L 7 + L 8 )(m 4 - 2m A K + m 2 K ml)) 

B.0.8 r/n — > r/ir 

This channel is pure 7=1 isospin. The amplitude is given by 

2 

— 771 

T 2 (s,t, U ): 



3fnfiT 



(54) 



T 4 (s,t,u) = 



{(t - 2ml)(t - 2m 2 r) )(6L 1 + L 3 ) + U L 4 (ml + 2m 2 K ) 



Vlf 2 

+ (3L 2 + L 3 

- 16m 2 K ml(L 4 - 2L 6 + 2L 7 ) - 3mlm 2 v (AL A + L 5 )} 



I 2 2\2 1 / 2 2\2 

(s -m n - m v ) + [u - m n - m^) 



+ m 4 (AL 4 -L 5 - 8L 6 + 32L 7 + 12L 8 ) 

(55) 
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Table I: Physical states used in the different /, J channels 





1=0 


1=1/2 


1=1 


1=3/2 


1=2 


J=0 


7T7T 

kk 


Ktt 
Krj 


711] 

KK 


Ktt 


7T7T 


J=l 


KK 


Ktt 
Krj 


7T7T 

KK 







Table II: Fit parameters Lj • 10 3 and comparison with the L\ ■ 10 3 of xPT 



q max = 1.02 GeV 




L 2 


4 




U 


2Lg ± Lg 


L 7 




0.5 


1.0 


-3.2 


-0.6 


1.7 


0.8 


0.2 


ji = 1.2 GeV 




rr 


^3 


r r 


rr 
^5 


2Lg ± Lg 






0.1 
±0.3 


0.9 
±0.3 


-3.5 
±1.1 


-0.7 

±0.5 


0.4 
±0.5 


0.0 
±0.3 


-0.4 
±0.2 
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Table III. Masses and partial widths in MeV 



Channel 


Resonance 


Mass 


Width 


Mass 


Width 


Partial 


(I, J) 


from pole 


from pole 


effective 


effective 


Widths 


(0,0) 


a 


442 


454 


« 600 


very large 


7171 - 100% 


(0,0) 


/o(980) 


994 


28 


« 980 


« 30 


7T7T — 65% 

KK - 35% 


(0,1) 


0(1020) 


980 





980 







(1/2,0) 




770 


500 


w 850 


very large 


Kn - 100% 


(1/2,1) 


K*(890) 


892 


42 


895 


42 


Ktx - 100% 


(1,0) 


a (980) 


1055 


42 


980 


40 


7T77 - 50% 

KK - 50% 


(1,1) 


p(770) 


759 


141 


771 


147 


7T7T - 100% 




800 800 

E cm [MeV] 





50 — 

: <5 n 



1100 

E cm [MeV] 



; '(d)' ' ' I 








\ 






















- (l-r?? )/4 







1100 

E em [MeV] 



50 1200 1000 1050 1100 1150 1200 

E cm [MeV] 

Fig.l: Results in the I = J = channel, (a) phase shifts for tttt — + irir as a fraction of the CM energy of the 
meson pair: full triangle (l9), open circle p0| , full square plj , open triangle j22], open square Q (all these 
are analysis of the same experiment [H ) , cross , full circle J25| , empty pentagon [^6| . (b) phase shifts 
for KK — > 7T7r: full square p7J , full triangle |2q]. (c) Phase shifts for KK — > i4T^T. (d) Inelasticity: results 



and data for (1 — 77 )/4: starred square |M, full square [E7| , full triangle |2S[|, full circle |29 
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Fig. 3: Phase shifts for tttt — > tttt in the I = 2, J = channel. Data: cross jjl], empty square |3^] 
triangle p3] ]. 
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Fig. 4: Results in the / — 1, J — channel, (a) phase shifts for ttt] — * tttj. (b) Invariant mass distribution 
for ttt] data from p3. (c) Phase shifts for KK — » ttti. (d) inelasticity. 
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Fig. 5: Results in the I = 1/2, J = channel, (a) phase shifts for Ktt — > if7r. Data: full circle J36]], cross 
p7| , open square full triangle j39fl , (b) phase shifts for Ktt — > if 77. (c) phase shifts for if 77 — > ifr/. (d) 
inelasticity. 
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Fig. 6: Results in the / = 1/2, J = 1 channel, (a) phase shifts for Kir — » Kir. Data: full triangle [ p6| , 
circle [pa. (b) phase shifts for Kit — > if 77. (c) phase shifts for Kr\ — > J^z/. (d) inelasticity. 
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Fig. 7: Phase shifts for Kir — ► i^7r in the / = 3/2, J = channel. Data: open triangle open circle 
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Fig. 8: (|T/j = oi|) 2 for A'if — > KK showing the singularity corresponding to the 4> resonance. 
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Fig 9: Imaginary part of the tttt amplitude in the (I,J)=(0,0) channel in the second Riemann sheet. On the 
left we show a three dimensional plot were we can observe the different structure of the a and /o poles. On 
the right we show a contour plot of the lower half plane of the second sheet. The a pole is very far away 
from the real (physical) axis and its lines of maximum gradient are parallel to it, in contrast with the /o- 
That is why the effect of both poles in the phase shifts (Figure 1) is so different. 




700 720 740 760 780 800 1000 1020 1040 1060 1080 1100 



Fig 10: The poles associated to the p (left) and ao (right) are oriented differently. The p mass seen on the 
(I,J)=(1,1) phase shifts is slightly bigger than the real part of the position of the p pole, whereas the peak 
of the mass distribution where the ao is observed (see Figure 4) is smaller than the real part of the ao pole. 
Concerning the widths, they are obtained as twice the imaginary part of the associated pole position. 
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